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We analytically treat the scattering of two counter-propagating photons on a two-level emitter 
embedded in an optical waveguide. We find that the non-linearity of the emitter can give rise to 
significant pulse-dependent directional correlations in the scattered photonic state, which could be 
quantified via a reduction in coincident clicks in a Hong-Ou-Mandel measurement setup, analogous 
to a linear beam splitter. Changes to the spectra and phase of the scattered photons, however, 
would lead to reduced interference with other photons when implemented in a larger optical circuit. 

We introduce suitable fidelity measures which account for these changes, and find that high values 
can still be achieved even when accounting for all properties of the scattered photonic state. 


The realization of on-chip all-optical information 
processing requires the implementation of quantum 
computation schemes in integrated photonic cir¬ 
cuits m, where the information is encoded in the 
travelling photons H E]. These schemes demand 
efficient single-photon sources mu, photon detec¬ 
tors la nm, and photon gates mum Since pho¬ 
tons do not inherently interact, in order to realize 
two-photon gates such as the controlled phase or the 
CNOT gate [lOj T3j[14], optical non-linearities or post 
selection schemes are required. As Kerr-type non- 
linearities are usually weak in the few-photon limit, 
promising candidates to mediate the two-photon in¬ 
teractions are quantum mechanical two-level-systems. 
When two photons scatter off a two-level system, non¬ 
trivial correlations can be induced in the scattered 
state [EJ!20j, the nature of which ultimately deter¬ 
mines the feasibility and scalability of optical circuits 
based on these components. Thus, a two-level-system 
embedded in a one-dimensional waveguide constitutes 
an important prototypical system in which to investi¬ 
gate few-photon scattering and non-linearity-induced 
photonic correlations. We note that such systems have 
been experimentally realized, for example, by self- 
assembled semiconductor quantum dots in photonic 
crystal waveguides, with emitter-waveguide coupling 
efficiencies reaching values in excess of 98% m- 

For two photons scattering on a single emitter, it 
is known that non-linearities are strongest when the 
photons are identical, and their spectral linewidths 
are comparable to that of the emitter, which results 
in the strongest correlations in the scattered photonic 
state mm- As such, identical input photons of a 
specific spectral lineshape are required. In the single 
photon case, however, it is know that such finite-width 
photons experience changes in both their spectra and 
phase as a result of the scattering process il 
which is also known to be the case for two-photon scat¬ 
tering m- Thus, once a photon has passed through 
an optical gate, it is no longer identical to an input 
photon, which in turn may limit the effectiveness of 
subsequent gates. This raises questions regarding the 
feasibility of integrating a large number of photonic 
gates needed to create complex optical circuits. The 


purpose of this work is to explore how two-photon 
pulses are altered by the scattering process, and inves¬ 
tigate how these alterations depend on the level of in¬ 
duced non-linearities. Interestingly, we find that non- 
linearaties can actually suppress spectral and phase 
changes, thereby increasing the similarity of the scat¬ 
tered and input photons. As such, even when correctly 
accounting for all properties of the scattered photonic 
state, fidelities between the scattered and a desired 
directionally entangled state as high as 80% can still 
be achieved. 

Various methods have been used to analyze multi¬ 
photon scattering in systems consisting of a local¬ 
ized scattering object coupled to a waveguide. These 
include fully numerical approaches da m as well 
as analytic approaches such as the input-output 
formalism ns, the real-space Bethe ansatz [24] . 
the Lehmann-Symanzik-Zimmermann formalism [25], 
Laplace transforms [26] . a wavefunction-based ap¬ 
proach nana, and master equation formalisms m- 
Several of these approaches allow for analytic deter¬ 
mination of single- and two-photon scattering matrix 
elements, which directly relate the scattered state to 
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FIG. 1: (a) Two counter-propagating single-photon pulses 
propagate toward a two-level system in its ground state, 
(b) The post-scattering state is measured by detectors in 
each chiral waveguide mode sub-group. 
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the initial state of the system. Studies of two ini¬ 
tially co-propagating photon pulses have been made 
for various scatterers coupling to waveguides, such as 
a single emitter mm, an emitter inside an optical 
cavitiy [25j EH] , and a non-linear optical cavity [29]. 

Here we give a largely analytical description of the 
scattering of two counter-propagating photons im¬ 
pinging on a two-level emitter in a one-dimensional 
waveguide, as sketched in Fig. [lja). We use the scat¬ 
tering matrix formalism, and analyse the strong non¬ 
linearity-induced correlations in the scattered state for 
various input states. In addition, we introduce fidelity 
measures to quantify the induced correlation in the 
scattered state, taking both the spectrum and phase 
of the scattered state into account, and discuss their 
experimental interpretations. 

This paper is organized as follows: In Section |T] we 
introduce our model. In Section|II|we review the scat¬ 
tering of a single-photon pulse on a two-level-emitter, 
and introduce fidelity measures to quantify the simi¬ 
larity between the incoming and scattered photons. In 
Section ED the formalism is extended to the scatter¬ 
ing of two counter-propagating single-photon pulses, 
where our fidelity measures are used to analyse in¬ 
duced correlations and spectral changes, and how 
these relate to the level of non-linearities. 


I. GENERAL THEORY 

We consider a quantum two-level-emitter coupled to 
two subsets of the modes in a one-dimensional waveg¬ 
uide as shown in Fig. [lj with the subsets differing by 
the direction of propagation. In the following we limit 
ourselves to lossless systems, though we note that this 
assumption could be relaxed by coupling our system 
to additional external reservoirs [33 33]. Additionally, 
we neglect waveguide dispersion in the considered fre¬ 
quency interval, and we assume a localized scatterer 
(dipole approximation), i.e. that the scattering occurs 
only at a single point in space. 

We describe the waveguide by two chiral mode sub¬ 
sets, being the right (mode index 1) and left (mode in¬ 
dex 2) propagating modes. The corresponding Hamil¬ 
tonian of the (bare) waveguide is 

2 />oo 

H 0 = V / dkhu(k)cl(k)ci(k), (1) 

i=l 

where each mode in subsystem 1 and 2 is characterized 
by a wavevector &, annihilation operator c*(fc) and en¬ 
ergy fou(k). In writing the Hamiltonian in this way we 
implicitly consider a single polarization of the waveg¬ 
uide modes. The Hamiltonian describing the emitter 
and its coupling to the waveguide is given by 

t 2 f‘Oo 

Hi = ~Y~cr z + hg^ / dk[a+Ci(k) + a-c\(k)} (2) 

where ujq is the resonance frequency of the emitter and 
g is its coupling strength to the waveguide, which is 


assumed to the frequency independent. This assump¬ 
tion is justified provided the linewidth of the emitter 
is small compared to the optical carrier frequencies of 
the photons. The operators cr + and <r_ are the cre¬ 
ation and annihilation operators of the emitter and 

(J z = (7_|_(7_ — <T_(T_|_. 

We consider pulses having the same carrier wavevec¬ 
tor k p and corresponding frequency u p = oj(k p ). It is 
therefore convenient to work in a frame rotating with 
this carrier frequency m- We relate the frequencies 
of the waveguide modes u(k) to their wavevectors us¬ 
ing a Taylor expansion around k p , giving 

u(k) uj p + v g (k - kp) (3) 

with v g = (du / dk)\ k=k being the group velocity. The 
rotating frame is defined by the transformation H 
H = UHW + i h(d t W)U with U = exp[iu p t(a z /2 + 
Ei /o°° dAct (jfe)ci(*))]. We find H = H 0 + Hi with 

2 poo 

H 0 = hv g Y, d kka\(k) 0 i(k), (4) 

where k = k — k p and we have defined the new anni¬ 
hilation operators ai{k) = c*(fc + k p ). The interaction 
Hamiltonian is now given by 

fe A 2 fOQ 

Hi = —a z + hg^2 / dk[a + ai(k) + a-a\(k)], (5) 

^ i= 1 ^ — 

with A = cjo — kpVg the detuning of the carrier fre¬ 
quency from the emitter transition frequency. We 
note that in obtaining Eqs. Q and (J5| we have ex¬ 
tended the lower limits of integration from — k p to 
—oo. This approximation is justified since we will be 
interested in pulses with wavevectors centered around 
k — k — kp = 0 and whose widths are much smaller 
than k p . 

II. SINGLE-PHOTON SCATTERING 

Before we consider the scattering of two photons, we 
first review the single-photon scattering case and in¬ 
troduce the scattering matrix formalism. We consider 
cases for which the emitter is initially in its ground 
state. Following Fan et al. [16], we relate the state 
of the system long after the scattering process to the 
initial state through the scattering matrix S. For non¬ 
linear scatterers the scattering matrix will in general 
be frequency dependent, and for a localized scatterer 
the scattering elements are defined as Pi 

i(p\S^\k)i = 2(p\S ( ' 1) \k) 2 = t k S{p-k), (6) 

2{p\S^\k)i = i(p\S^\k) 2 = f k 5(p-k), (7) 

where the notation implies, a\(k)\(f>) = \k)i with \<f>) 
the vacuum, and where tk and fk are the frequency- 
dependent single-photon transmission and reflection 
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coefficients, respectively. The delta-functions reflect 
momentum conservation, and as no external loss chan¬ 
nels are present, \tk\ 2 + \ fk\ 2 = 1. 

An arbitrary single photon state propagating to the 
right is written 

/ oo 

dkak)a\(kM, (8) 

-oo 

where £(fc) is the wavepacket in momentum space. 
The post-scattering state corresponding to the incom¬ 
ing state expressed in Eq. (J8| is defined as |£)t_>.oo = 
|<$;o)- It is obtained by inserting the identity oper¬ 
ator 

/ oo 

dp\p)ii{p\, ( 9 ) 

-°° 

from which we find 

/ OO /‘OO 

dpt P t(p)\p)i+ dpf p ^(p)\p) 2 . 

-oo J — oo 

( 1 ") 

The two terms above reflect the fact that the 

photon can be transmitted or reflected. The 

scattering probabilities are defined as Pi = 

t—YOO (£| X-^oo |p)i»(p|£)t— k»> With the transmission 
and reflection probabilities corresponding to i = 1 and 
% — 2 respectively. These probabilities are found to be, 

/ OO /-OO 

dp\ipt(p)\ 2 , P2= d P\fpti(p)\ 2 - 

-oo J —oo 

(n) 

The theory above applies to any localized scatterer 
interacting with two chiral waveguide modes. We 
now specifically consider the emitter-waveguide sys¬ 
tem sketched in Fig. [T] and described by the Hamilto¬ 
nians in Eqs. @ and (j5]). In this system the reflection 
and transmission coefficients may be found through 
calculation of the single-photon scattering matrix ele¬ 
ments m, which gives 

= k -A = ~if/ (2vg) 

k k — A + iT/(2vg) ’ k k — A + \T/{2v g ) ’ 

( 12 ) 

where T = 47 rg 2 /v g is the decay rate of the emitter. 
We note that this form of T can be found through 
Fermi’s Golden Rule, and is valid for a lossless sys¬ 
tem in which the emitter couples equally to modes 
propagating in both directions. 

We consider three transform-limited incoming pulse 
shapes, i.e. Lorentzian, Gaussian, and step- 
function OHS, which are defined by the wavepackets 


^Lor(^) — 

vV( 27 0 

k — i<r/2 ’ 

(13) 

CGauss(^) = 


(14) 

Csquare(^) 

a~V 2 e{a/2- |fc|), 

(15) 


with a' = (2yffii(2)) -1 a for the Gaussian wavepacket, 
and where 0(k) is the Heaviside step func¬ 
tion. All wave packets are normalised such that 



z in units of vJT 

FIG. 2: (a) Intensity spectra, |£(fc)| 2 , for the three single¬ 

photon wavepackets in Eqs. |l3])-([l5]), here plotted with 
a = r /v g . (b) Corresponding spatial pulse profiles, |£(z)| 2 , 
with large values of z corresponding to the front part of 
the pulse, arriving first at the position of the emitter. 


dfc l£(^)| 2 = 1? and have a spectral full width- 
half maximum of a. The intensity spectra of the 
pulses are shown in Fig. [2] together with the spa¬ 
tial pulse profiles, defined as the Fourier transform 
t(z) = (27r)~ 1 / 2 /^ o dfc^(fc)exp[ifcz]. 

A comparison of the resulting reflection probabili¬ 
ties is shown in Fig. [3ja) for each of the three single¬ 
photon wavepackets in Eqs. (p~3]) - ([15] ) , as has been 
calculated in previous works [171122) . The frequency 
components of the pulse closest to the transition en¬ 
ergy of the emitter interact most strongly, and those 
at the exact frequency of the emitter (k = 0) are 
perfectly reflected [22]. Thus, as the spectral pulse 
width is decreased, the reflection probability increases, 
and reaches unity for resonant monochromatic pulses 
(a —> 0). In the opposite limit of a —> 00 , only a 
vanishing fraction of frequency components overlap 
with the spectrum of the emitter, resulting in com¬ 
plete transmission since the pulse does not interact 
with the emitter. The pulse shape also has an impor¬ 
tant impact on the reflection probability. Since the 
Lorentzian has the largest spread of frequency com¬ 
ponents for a given FWHM, it interacts least with the 
emitter and correspondingly results in the lowest re¬ 
flection probability. 

If many emitters are to be implemented in a larger 
sequence of photonic devices or gates, it is important 
that scattered photons maintain their spectral prop¬ 
erties, i.e. the pulse shape and phase variation across 
the pulse. We therefore seek to define measures to 
compare the scattered state with some desired out¬ 
put state. We define three measures for this purpose, 
each with a different physical significance. We use the 
quantum state fidelity m as a measure of the degree 
to which the scattered state is quantum mechanically 
identical to the desired state (neglecting overall phase 
differences), 


*full — | (£des |£)t— >00 | 


(16) 
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FIG. 3: (a) Reflection probabilities as a function of pulse 

width for the three pulse shapes given in Eqs. (fl3]) - (fl5]) 
with carrier frequencies resonant with the emitter, (b) The 
three fidelity measures from Eqs. (17)-(21) for scattering 
of a Lorentzian pulse shape. 


where |£des) is the desired state. When the desired 
state is equal to the input state but propagating in 
the opposite direction (i.e. a reflected but otherwise 
unchanged state, |£ des ) = /^ c C(p)b) 2 ), we find 


Fm 



dpr(p)\£(p)\ 2 


2 


(17) 


caused by the absorption in the emitter. Using the 
defined Fourier transform and Holder’s inequality, we 
find Fi n t is an upper bound for this spatial fidelity, 


-^spat — max 

oz 


r 


dp\t(p)\ 2 r(p)e' 


i8zk 


2 

< Fi nt . 


( 20 ) 


Finally, when neither the phase nor the spectral dis¬ 
tribution are important, the scattered state may be 
projected onto a basis which merely counts the num¬ 
ber of photons in each waveguide mode, e.g. as in 
Fig. 0b). The fidelity in this case simply becomes 
the probability of detecting a photon in the desired 
output mode (here the reflected field), 

^prob = Pr • (21) 


We evaluate these fidelities for a Lorentzian input, 
and show the results in Fig. [3|b). For Ff u n and F pro b 
we find the exact expressions 


*full 

Tprob 


CO 2 

(f + a) 2 +4A 2 ’ 

(f+g)f 

(f + a) 2 +4A 2 ’ 


( 22 ) 

(23) 


where f = T/v g . From Fig. [3^b) we see that 0 < 
Ff u ii < Fi n t < Fp r ob < 1. This reflects the progres¬ 
sively less stringent criteria of these three measures. 
As the desired state in each case is a fully reflected 
state, the fidelities are largest for small FWHMs. 


III. TWO-PHOTON SCATTERING 


This fidelity measure is appropriate, for example, if 
the scattered photon were to interfere with a second 
input photon in a Hong-Ou-Mandel-type interference 
experiment (where two truly indistinguishable pho¬ 
tons exit a 50/50 beam-splitter in the same arm). 

For cases in which the phase of the scattered state 
is unimportant, but instead we are interested in how 
the intensity in the pulse is distributed spectrally, the 
similarity between the scattered and the desired pulse 
may be characterized by 

^int= (/ , ^dp|f(p)||e(p)| 2 ) 2 . (18) 


This fidelity measure would be relevant when com¬ 
paring the energy distributions in the scattered and 
desired pulse, which could be achieved by introducing 
spectrometers in a setup as sketched in Fig. 0b), but 
disregarding the arrival times at the detectors. It can 
be seen to be the limiting form of a spatial definition 
of a fidelity measure, defined as 


T sp at = max 

oz 


f 


dzC(z)€acat(z ~ Sz) 


(19) 


where £ S cat(^) is the spatial representation of the scat¬ 
tered wavepacket, which after the scattering may be 
displaced by Sz in the rotating frame due to a delay 


We now turn to the main focus of this work, and 
extend our formalism to describe the scattering of two- 
photon states. In the single-photon case, energy con¬ 
servation implied that an approximately monochro¬ 
matic single-photon wavepacket would scatter without 
changing its frequency. In the two-photon case, energy 
conservation only demands that the sum of the ener¬ 
gies of the two incoming and two scattered photons 
is conserved. According to Fan et al. m, we can de¬ 
fine a two-photon scattering matrix, S^ 2 \ in a similar 
way to S^\ which contains terms describing single¬ 
photon scattering, and also additional terms stem¬ 
ming from two-photon scattering processes. The ad¬ 
ditional terms involve four-wave mixing mechanisms 
between the two incoming and two scattered pho¬ 
tons Usuis]. 

In the rotating frame, a general two-photon state in 
the momentum representation is written 


II 

C. 

POO 

/ d k 

' — oo 

POO 

/ dk! Pu(k,k/)al(k)a\(kf)\tf>) 

' —oo 

+ 7I u 

POO 

/ dA; 

' — oo 

POO 

/ dk'/3 2 2(k, k')a\(k)a\(k')\<f>) 

1 — OO 

A 

POO 

/ d/c 

' —oo 

POO 

/ dk' f3rz(k, k')a\(k)a\(k')\<f>) 

! — oo 


( 24 ) 
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normalized such that d k d k! (|/3n(fc, k')\ 2 + 

|/? 12 (k,k')\ 2 + \(322(k,k')\ 2 ) = 1. Introducing the no¬ 
tation | kk!)a> — {i'i(k'k\y = aj(fc)at,(fc')|0), the two- 
photon scattering elements are ns 

J y(pp'|S' (2) |fcfc') i ^ = cxji } kOij>i> ! k'S(k — p)S(k' — p') 

+ oty i}k otji' }k ’5(k - p')S(k' - p ) 

+ j Bpp>kk'S(p + p' -k- k'), 

(25) 

for i,j G {1,2}, and where 

4 if i = j 


OLj — 


r k if i ± j 


(26) 


are the single photon reflection and transmission ma¬ 
trix elements. Here B pp >kk' describes interactions be¬ 
tween the two incoming and the two scattered pho¬ 
tons, and is determined by the specific localized scat- 
terer considered. As in the single-photon case, to find 
the scattered state, we insert the identity operator, 
which is now given by 


U 


/ OO POO 

dp / dp' 

-oo J — OO 


ib'p)llll(pp'| + |p'p>2112(pp'| + ^\p'p)2222{pp'\ 

(27) 

If we assume an initial state consisting of two counter- 
propagating photons, (3 12 is the only non-zero expan¬ 
sion coefficient in Eq. (24), which results in the post¬ 
scattering state 


/£—>- oo 


/ oo poo r 

dp / dp' 

-oo J —oo l 


itpf p f +f p i p ')Pi 2 (p,p') + \b 12 {p,p') 
(‘ r p r p > +i p ip')p 12 (p,p') + ^b 12 (p,p') 
(t p f p , + r p t pl )p 12 (p,p') + \b 12 (p,p') 


a\{p)a\(p') 

a i(p) a Up') 

4(p)4(p') 

<£)• (28) 


The first term in each of the three square brackets 
in Eq. (28) represents single photon scattering pro¬ 


cesses, where the first factor contains the appropri¬ 
ate combinations of transmission and reflection coef¬ 
ficients connecting the initial and final photon con¬ 
figurations. Multi-photon process are contained in 
the pulse-dependent contribution which de¬ 

scribes processes induced by the emitter non-linearity, 
and is given by 

/ oo 

dk/3 12 (k,p+p' - k)B pp , k(p+p ,_ k) . 

-OO 

(29) 


We define Pu (P 22 ) as the probability that both 
photons are measured propagating in waveguide mode 
1 (mode 2), and P 12 the probability that one photon 
propagates in each waveguide mode. From Eq. (28) 
we find 


Ai = - 


1 p 00 p 00 

- d P dp' 

^ J—00 J—00 


( t p r p / +r p t p ’)(3 12 (p,p') + -bi 2 (p,p') 


Pm = 


/ OO POO 

dp / dp' 

-00 J — 00 


C Wp' + r p r p >)Pi 2 (p,p') + ~.b 12 {p,p') 


(30) 


(31) 


with P n = P 22 and Pu + P 12 + P 22 = 1. 

As in the single photon case, we will be interested in 
comparing the scattered two-photon state described 
by Eq. (28) to some desired state using the fidelity 
measures we have introduced in Eqs. ([T7|)-(|2l|). For 
this purpose we consider the scattered two-photon 
state that would be obtained if the scatterer were re¬ 
placed by a perfect 50/50 beam-splitter, which pre¬ 
serves the shape and phase of the input photons. In 
this case the scattered state would be 

1 POO POO 

|/?des} = ^/ d k dk'P 12 (k,k') 

V ^ J — oo J — oo 


—=a\(k)a\(k') + —=al(k)al(k ') 


). (32) 


L v/2 1W 1V ' y/2 

With this desired state our fidelity measures become 

2 


*fdll = 


/ OO POO 

dp / dp'/3{ 2 (p, p')/3{ 2 (p, p') 

-oo J — oo 


1 —oo J —oo 

POO POO 


/ oo poo i 2 

dp dp'|/3i 2 (p,p')| \Pi 2 {piP')\ 

-oo J— oo 


Pint — 

LJ — oo J—oo 

Ppi-ob = Pll + P 


22 


(33) 

(34) 

(35) 


with /3[ 2 (p,p') = (t p r p i + r p t p >)P 12 (p,p') + \b 12 {p,p'). 


A. Two-level scatterer 

The theory presented above is valid for any localized 
scatterer. We now specifically consider the two-level- 
emitter-waveguide system described by Eq. ([5]), and 
focus here only on pulses starting equidistantly from 
the emitter. Furthermore, we only treat pairs of in¬ 
put pulses with the same spectral linewidth, although 
the formalism can be straightforwardly extended to 
more general cases. For a two-level-emitter the single 
photon transmission and reflection matrix elements 
tk and ffc are given by Eq. (12), while the two-photon 
scattering element is [16] 


d .>/r 

Ppp'kk' — i SpSp'y^k + 

7r 


(36) 
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where 


$k — 


Vf/i 


k - A + ir/(2i; s )' 


(37) 


We only consider uncorrelated photon input states, 
and as such /3(k,k') is a symmetrized prod¬ 
uct of two single-photon wavepackets /3(k,k') = 

+ £'(&)£(&')], which is normalized as 

r o od*/“cod* , l/3(*,* , )l a = l- 

We begin our analysis of the scattered state by con¬ 
sidering correlations in photon detection events in the 
two waveguide mode subsets, as depicted in Fig.[ljb). 
In this case no information regarding the spectrum 
and phase of the scattered photons is obtained, and 
the appropriate fidelity measure is F pro b, which is 
equal to 1 minus the probability of detecting a co¬ 
incidence in the two detectors, i.e. for F pro b = 1 no 
coincidence events are measured (a perfect Hong-Ou- 
Mandel dip would be observed). In Fig. |4^a), F pro b is 
calculated for Gaussian and Lorentzian input pulses 
for zero detuning (A = 0), for which the interaction 
between the pulses and emitter is greatest. We see 
that very high fidelities are obtained, reaching values 
of ^ 80% for the Lorentzian input and ~ 90% for the 
Gaussian. Maximal correlations are achieved in the 
regime where the emitter and pulse linewidth are sim¬ 
ilar, as has been demonstrated numerically in earlier 
work m- Interestingly, although the Lorentzian pulse 
shape is well-known to be the optimal pulse shape for 
maximally exciting the two-level emitter with a single 
photon [35!, it is not the optimal shape for maximizing 
the directional correlations in the scattered state. 

The high fidelities obtained demonstrate that the 
scattered states are highly directionally entangled, in 
analogy with the effect of an optical beam-splitter. 
However, in contrast to the classical beam splitter, 
the high correlation seen here is induced solely by 
non-linearities. To demonstrate that the high direc¬ 
tional correlations indeed stem from non-linearities, 
in Fig. ga) we also show the case where the non¬ 
linear two-photon interaction term 612 (p,p') has been 
artificially set to 0 (dashed curves). For an uncorre¬ 
lated two-photon input pulse which is resonant with 
the emitter and which has a symmetric spectral wave- 
function amplitude, |£(— k)\ 2 = |£(£;)| 2 , as is the case 
here, Eq. (30) reduces to 


Pn = a(l - a), 


with a = 


r 


dp 


(f/2) 2 


p 2 + (r/2) 2 


\m\ 2 , (38) 


showing that Pn maximally attains the value 1/4, oc¬ 
curs when a = 1 / 2 . Thus, P pro b = 2 Pn never exceeds 
1 / 2 , as confirmed in Fig.[4^a), which indicates that no 
directional entanglement is present in the scattered 
state m- For 5 i 2 (p,p / ) = 0, the emitter behaves 
as a linear component (e.g. a lossless optical cavity) 
and cannot mediate interactions between the two pho¬ 
tons. As such, the scattering process is determined 



FIG. 4: Degree of directional entanglement P pro b 

plotted for varying incoming pulse widths, shown for 
both Lorentzian and Gaussian input pulses (solid lines). 
The corresponding values obtained when neglecting two- 
photon scattering terms are also shown (dashed lines). (a) 
Pulse and emitter are resonant, A = 0. (b) Pulse and 

emitter detuned by A = F/(2v g ). 


entirely by interference effects, which, unlike an opti¬ 
cal beam-splitter, cannot create entanglement in this 
system when the pulses are resonant with the emitter. 

A more direct analogy with a 50/50 beam split¬ 
ter can be obtained by detuning the input pulses by 
half the emitter linewidth, A = T/( 2v g ). For this 
value of the detuning, a monochromatic single-photon 
pulse will be reflected/transmitted with 50% prob¬ 
ability (whereas at A = 0 a monochromatic single¬ 
photon pulse is fully reflected). Fig. [4^b) shows P pro b 
for A = T/( 2 v g ), and we confirm that P pro b 1 as 
a —)> 0 as expected. The change in the fidelity due to 
the non-linearities now becomes smaller than in the 
resonant case, as the interaction between the pulse 
and the emitter is less efficient off resonance. Inter¬ 
estingly, in this case, for small <7, the non-linear inter¬ 
action actually deteriorates the beam splitting effect, 
since now the directional entanglement can be gener¬ 
ated by interference effects only. 

For a Lorentzian input analytic expressions for 
PiqLor = P 22 , Lor may be derived. We find 


„ 3fa(3a+f)(a+f)+4A 2 f(a+2f) , on , 

-Pll, Lor = T-1-=TT-1-T 


(3<j + f ) 2 + 4A 2 (a + f ) 2 + 4A 2 


with f = T/vg. The maximum value on resonance 
(A = 0) is obtained for <r/f = 3 -1 / 2 ~ 0.57, at which 
point -PipiLor ^ 0 *d 0 (and P P rob — 2 Pii 5 Lor ~ 0 . 8 ), in 
agreement with Fig.pFa). In comparison, with no non¬ 
linear terms, 612 (p,P) = 0 , the scattering probability 
becomes 


_ <7f(<7 + f) 2 +4A 2 f(<7 + 2f) 

“ 11 , Lor — r I 72 


(<r + r ) 2 + 4 A 2 
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p [units of T/vJ 


p [units of T/vJ 


p [units of T/vJ 



p [units of T/vJ p [units of I 7 v g ] p [units of T/vJ 


FIG. 5: Intensity spectrum of an incoming Gaussian two-photon state for two counter-propagating photons (left column), 
and the resulting two-photon intensity spectra for the scattered state (middle and right columns) for photons scattering 
in different directions, 12, and where both photons propagate in the same direction, 11 (identical to 22), with A = 0. 
The spectral width of the input pulses is varied: a = 0.2 T (1st row), a — IT (2nd row), and a — AT (3rd row). The 
intensity spectra are also shown for scattering with the non-linearity turned off, 612 (p, p) = 0, using a = 1 T (4th row). 


which has an on-resonance maximum for a = T, giv¬ 
ing Fp r ob = 2 Prrlqt = 1/2. Thus, for the Lorentzian 
pulse, non-linearities increase the scattering probabil- 
ity by a factor of A 1 , Lor /^T, Lor = 1 + 2 /(l + 3<r/f) 
on resonance. For a -T oc, the interaction with the 
emitter becomes infinitely weak and no enhancement 
is present. In the opposite limit of a 0, the en¬ 
hancement factor is 3. 


B. Scattering Fidelities 

As discussed above, if the scatterer is to be im¬ 
plemented in a larger optical circuit, in addition to 


considering in which direction the photons scatter, 
the amplitude and phase of the different frequency 
components may also be important. In such a case, 
Fp r ob is no longer a sufficient fidelity measure, since 
it contains only directional information. As illus¬ 
trated in Fig. [5j where we plot the intensity spec¬ 
trum of a scattered Gaussian wavepacket, the spectra 
of the scattered pulses change significantly during the 
scattering process. For input pulses with a narrow 
spectral linewidth compared to the emitter (1st row), 
the pulse power at the emitter position remains low 
due to the corresponding broad spatial profiles of the 
pulses. In that case, the non-linearity is only weakly 
addressed, and the individual photons are predomi- 





























1 



FIG. 6: Fidelities from Eqs. (33)-(35) plotted for varying 
width of the input pulses, shown both for Gaussian and 
Lorentzian inputs on resonance with the emitter. 


nately reflected. A weak non-linearity-induced four- 
wave mixing process is signified by the appearance of 
diagonal features, as one photon achieves a larger en¬ 
ergy, while the energy of the other decreases. When 
the pulse and emitter linewidths are comparable ( 2 nd 
row), the predicted strong directional correlation is 
induced HE with the pulse profile being almost pre¬ 
served. For spectrally broad pulses (3rd row), only 
the near-resonant part of the spectrum interacts with 
the emitter. We see that the spectral components at 
the emitter frequency are absent from the transmitted 
pulse since these have been reflected without signifi¬ 
cant two-photon effects. 

Interestingly, the fact that the pulse spectrum is 
almost perfectly preserved when the pulse and emit¬ 
ter linewidths are comparable ( 2 nd row) can be at¬ 
tributed to non-linearities. This can be seen in the 
4th row, where we again show the initial and scattered 
spectra for the case a = T, but where we have artifi¬ 
cially set the non-linear term equal to zero, 612 = 0 . 
By comparison with the 2nd row, we can clearly see 
that the non-linearities not only give rise to the direc¬ 
tional entanglement, but also suppress changes to the 
spectral shape. 

To quantify both the spectral and phase deviations 
between the scattered and the desired state, all three 
fidelities defined in Eqs. (33)-(35) are shown in Fig. [ 6 j 
By comparing F pro b to Fi nt , i.e. taking into account 
the difference in the spectra of the scattered and de¬ 
sired pulse (but not the phase), we see that the fi¬ 
delity becomes lower, and most significantly so for 
pulses with a small spectral linewidth. This can be 
understood from Fig. [5j where we see that the scat¬ 
tered wavepacket for the spectrally narrow input ( 1 st 
row) is clearly influenced by strong four-wave mixing 
effects. For pulses with larger widths, these effects 
are weaker, since a larger fraction of frequency com¬ 
ponents are detuned from the emitter transition and 
therefore interact only weakly. 

Scattering-induced phase differences across the 
pulses may be examined by comparison of F[ nt and 
Ff u n. As illustrated in Fig. [ 6 j these fidelities are al¬ 


most equal for spectrally narrow pulses, whereas sig¬ 
nificant deviations are seen for spectrally broad pulses. 
This may be explained by considering the simpler 
single-photon scattering case. From Eq. (Jl2|, we see 
that a resonant, monochromatic pulse will be reflected 
with a phase shift of 7 r, whereas spectrally broader 
pulses attain a phase shift from n/2 to 3tt/2 across 
the pulse spectrum. Thus, spectrally broad pulses ex¬ 
perience larger decreases in the fidelity due to phase 
mismatching with our given desired state. 

As the phase changes correspond to modifications 
to the spatial profile of the pulse, we specifically con¬ 
sider how the spatial pulse profile (here analogous to 
the temporal shape) is changed during the scattering 
process. To clearly illustrate the effect of the non¬ 
linear scattering on the spatial profile, we evaluate 
the photon density at a specific point of the photon 
wavepacket in the rotating frame, defined as 

N(z) = t ->00 (/3|a f (z)a(z)\/3) t ^oo (41) 

where a(z) = (27t) -1 / 2 f^° dfca(fc) exp[ifc^] is the an- 
nihilation operator for an excitation at a position z in 
a frame rotating with the pulses. For Gaussian and 
Lorentzian input pulses, the photon density is plot¬ 
ted in Fig. [7| For both pulse shapes, we see that a 
delay occurs due to interaction with the emitter, and 
furthermore the non-linearity improves the similarity 



z [units of Vg/T] 



FIG. 7: The photon density in the moving frame, N(z), 
calculated for the input pulse and for photons in the same 
(11 identical to 22) or in different modes after the scat¬ 
tering (12) for resonant input pulses with a — IT and 
A = 0. The largest values of z correspond to the front part 
of the pulse, and the solid (dashed) lines include (do not 
include) the non-linear two-photon interaction, bi 2 (p,p). 
(a) Gaussian wavepacket (b) Lorentzian wavepacket. 
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between the scattered and incoming field. The Gaus¬ 
sian pulse is observed to preserve its spatial symme¬ 
try, as compared to the Lorentzian input pulse. This 
is due to the fact that the part of the photon pulses 
which is absorbed by the emitter is re-emitted with an 
exponential shape that is spatially reversed compared 
to the input pulse, which explains why F pro b deviates 
significantly from F[ nt for large spectral linewidths in 
Fig. |6| 

IV. CONCLUSION 

We have analytically demonstrated that the non¬ 
linearity of a two-level-emitter can induce strong 
pulse-dependent directional correlations (entangle¬ 
ment) in the scattered state of two initially counter 
propagating photons. These correlations are maxi¬ 
mized for photons with spectral widths comparable 
to that of the emitter, and also depend on the spe¬ 
cific spectral shape of the photons. Furthermore, we 
have investigated how the spectra and phase of the 
photons are affected by the scattering process, and 
introduced different fidelity measures to quantify the 
similarity of the scattered and input photons. Inter¬ 
estingly, for photons with spectral widths comparable 


to the emitter linewidth, where the directional corre¬ 
lations are maximized, the non-linearity of the emitter 
acts to suppress changes in the spectra and phase of 
the photons. As such, even when taking all proper¬ 
ties of the scattered state into account, a comparison 
to perfect directionally entangled photons with pre¬ 
served spectra and phases gives fidelities as high as 
~ 80% for Gaussian pulse shapes. A comparison of 
our fidelity measures indicates that when engineering 
photonic gate structures and other functionalities us¬ 
ing two-level-emitters, it is important to also consider 
spectral and phase changes when determining the ef¬ 
ficiency and scalability of non-linear photonic devices. 
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